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Thermal noise of a mirror is one of the limiting noise sources in the high precision measurement such as gravitational-wave detection, 
and the modeling of thermal noise has been developed and refined over a decade. In this paper, we present a derivation of coating 
thermal noise of a finite-size cylindrical mirror based on the fluctuation-dissipation theorem. The result agrees to a previous result with 
an infinite-size mirror in the limit of large thickness, and also agrees to an independent result based on the mode expansion with a thin- 
mirror approximation. Our study will play an important role not only to accurately estimate the thermal-noise level of gravitational-wave 
detectors but also to help analyzing thermal noise in quantum-measurement experiments with lighter mirrors. 
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1 Introduction 

Advancement in the reduction of technical noise and isola- 
tion of seismic vibration has let a high-precision measure- 
ment device like the interferometric gravitational-wave 
detectors [1] be so sensitive that tiny thermal fluctuation 
of the measurement surface can limit the sensitivity. It is 
important to develop a method to estimate the thermal- 
noise level as accurately as possible. Our study with a 
finite-size mirror is an upgrade of previous works for coat- 
ing thermal noise with some approximations. For mir- 
rors currently planned to be used in the next-generation 
gravitational-wave detectors, the difference is a few per- 
cent between the results with our finite-size analysis and 
with a conventional infinite-size analysis. Besides, our 
analysis, for the first time, includes the effect of tempera- 
ture fluctuations in the substrate and the coatings coher- 
ently summed up, with which the estimation of thermal 
noise will be more accurate at lower frequencies. The 
difference from the previous results can be larger if the 
mirror is thin. Thermal noise is also important in cold 
damping experiments [2] , where the purpose is to reach a 
quantum limit with a low-mass mirror, which may tend 
to be thin. In this paper, we show calculation results of 
coating thermal noise with a broad range of aspect ratios, 
which agree to the previous results with an infinite-size 
mirror in the thick limit, and also agree to the results 
with a thin-plate that are calculated using the modal- 
expansion method. It is important to know the thermal- 
noise level in the middle range so that an appropriate 
mirror can be used in the various experiments. 

There are two different ways that non-zero tempera- 
ture causes fluctuation of the surface of a mirror. The first 



one is via volume fluctuation under fixed temperature; 
called Brownian thermal noise. Brownian thermal noise 
in the power spectrum density (m/vTIz) is proportional 
to the square-root of temperature VT, besides the me- 
chanical loss has some temperature dependence. The sec- 
ond one is via temperature fluctuation that converts into 
the surface fluctuation through the thermal expansion 
and through the change of the refraction index. Ther- 
mal noise through the expansion is called thermoelastic 
noise [3] and thermal noise through the change of the re- 
fraction index is called thermorefractive noise [1][S]; the 
coherent sum of thermoelastic noise and thermorefractive 
noise is called thermo-optic noise in Ref. [BJ. Thermo- 
optic noise in the power spectrum density is linearly pro- 
portional to T, besides some parameters like the thermal 
conductivity or the thermal expansion depend on the tem- 
perature. 

Brownian thermal noise is related to the mechanical 
loss angle <f>. A current gravitational-wave detector em- 
ploys a mirror made of silica coated by tantala-silica dou- 
blets, and the loss angle of the silica substrate is several 
orders lower than that of the coatings [7] [S]. Thermo- 
optic noise is related to the heat flow in the r-direction 
(transverse to the beam) and in the z-direction (along 
the beam) of the cylindrical mirror. The both contri- 
butions are to be taken into account in the case with- 
out coatings [9], while the latter becomes dominant with 
coatings according to the difference of the mechanical pa- 
rameters of the materials |10j . In this paper, we focus on 
the derivation of Brownian thermal noise in the coatings 
and thermo-optic noise in the z-direction. 

Historically, Brownian thermal noise of a mirror had 
been analyzed using a so-called modal-expansion method . 



Gillespie and Raab demonstrated a calculation with Hutchintic equation, and the heat equation. In Sec. El we use 



son's method to derive the contribution from each me- 
chanical mode of an axisymmetric cylinder [12] . The 
contributions are added up with a weight function given 
by the power distribution of the Gaussian beam that 
probes the mirror. In 1998, Levin proposed a new way to 
analyze thermal noise using fluctuation-dissipation the- 
orem [13] I14j . Thermal noise is given from the multi- 
ple of the loss angle and the elastic energy of a mirror 
imaginarily pushed by the Gaussian beam. Levin demon- 
strated a calculation for substrate thermal noise with an 
approximation that the mirror is an infinite half-space, 
which is reasonable as the beam size on the mirror is 
usually set sufficiently small compared with the mirror 
radius to avoid the diffractive loss. Bondu et al calcu- 
lated substrate thermal noise of a finite cylinder using 
Levin's method [15]. The elastic energy of a monolithic 
substrate was derived with the boundary conditions of a 
finite cylinder. Harry et al extended the elastic equation 
with coating layers on the substrate and calculated Brow- 
nian thermal noise of thin coatings on an infinite mirror. 
In this paper, we will derive Brownian thermal noise of 
coatings on a finite cylinder. 

Thermoelastic noise, introduced in Ref. [H] , is related 
to the thermal expansion, and the fluctuation-dissipation 
theorem is again useful to derive the expansion. Liu 
and Thorne calculated thermoelastic noise of an uncoated 
substrate associated with the heat flow both in the r- 
and z- directions using Levin's method [9]. Thermoe- 
lastic noise of a coated material by the heat flow in the 
^-direction was calculated by Braginsky et al with an ap- 
proximation that thermoelastic dissipation via non-zero 
relaxation time of the heat flow at the coatings be re- 
garded instantaneous - thin- coating approximation [3]. 
They calculated thermoelastic noise both in the infinite 
case and in the finite case with the thin-coating approx- 
imation. Fejer et al calculated thermoelastic noise with- 
out the thin-coating approximation; the heat equation 
was solved both in the substrate and in the coatings, but 
the mirror was an infinite half-space and the expansion 
was approximated to be constant in z 10J. In this paper, 
we will derive thermoelastic noise by the heat flow in the 
z-direction without the thin-coating approximation and 
with a finite cylinder. 

Thermorefractive noise is calculated, even in this pa- 
per, with the thin-coating approximation. More rigorous 
analysis could be done but it will require an individual 
treatment of each layer with as many boundary condi- 
tions as the number of layers (typically 15 ~ 40), which 
shall be remained as a future work. Nevertheless, as will 
be shown in this paper, thermorefractive noise is not as 
sensitive to the thickness of the substrate as other two 
noise sources. 

The structure of this paper is as follows. In Sec. 
we explain the fluctuation-dissipation theorem, the elas- 



Bondu's solution to the clastic equation with a finite cylin- 
der and extend it with the coatings to calculate Brownian 
thermal noise of the coatings. In Sec.[H we extract the ex- 
pansion term from the last result and put it into the heat 
equation to calculate thermoelastic noise. In Sec. El we 
show the heat equation for thermorefractive noise, and 
the result will be combined with thermoelastic noise to 
make thermo-optic noise, which is shown in Sec. El In 
Appendix [A] we show the results for Brownian thermal 
noise and thermoelastic noise calculated with the thin- 
plate approximation, which should agree to the results of 
the finite-mirror calculations in the thin limit. In the end, 
Appendix |B] is a list of the parameters. 

2 Overview of the method 

2.1 fluctuation-dissipation theorem 

The conventional modal-expansion method and Levin's 
method are substantially equivalent methods to derive 
the noise spectrum using the fluctuation-dissipation the- 
orem. In both methods, an imaginary force is applied 
to the mirror. While the modal-expansion method first 
calculates the thermal motion of the mirror in many elas- 
tic eigenmodes and then adds them up with a weighting 
function for the Gaussian beam, Levin's method directly 
calculates the dissipation and thermal noise without the 
modal decomposition. In the main body of this paper, 
we use Levin's method. We also use the modal-expansion 
method in Appendix. [A] to calculate thermal noise of a 
thin plate. Since these two methods are quite different, 
coincidence of the results in the thin limit validates our 
calculation. 

The equality of the fluctuation and the dissipation is 
the important part of the fluctuation-dissipation theorem; 
the power spectrum of thermal motion is expressed by the 
following equation: 
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where fee is the Boltzmann constant and Z(Cl) is the 
impedance of the system, which is given by 
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The phase difference <p between the imposed force F = 
F cos (fit) and the resulting motion x = xo cos (fM — (f>) 
is called loss angle. The averaged dissipated power is the 
product of F and x in the same phase: W — FoXo£l(f)/2, 
so, with Eq. the power spectrum is rewritten as 



s x (n) = 



8k B TW 



8fc B T 



U<f> 



(3) 



2 



Here U is the maximum elastic energy that can be gener- 
ated by the imaginary force. The dissipation for Brown- 
ian thermal noise is derived with the elastic equation, and 
the dissipation for thermoelastic noise is derived with the 
heat equation. 

The logic above should be retraced in the case of ther- 
morefractive noise, although Eq. ([3|) still works as well. 
Thermorefractive noise is no actual motion but phase 
shift of the light due to the change of refraction index. 
Thus, it is not the imaginary force but imaginary entropy 
that should be added to the system to use the fluctuation- 
dissipation theorem. Adding the imaginary entropy, we 
can derive the noise spectrum of the temperature fluctu- 
ation, which can be converted into the fluctuation of the 
optical length. See Ref. [14] for the detail. Note that 
we regard the imaginary force, or imaginary entropy, to 
be static. This should be good if the frequency of our 
interest is much lower than the resonances of a mirror. 

2.2 Elastic equation 

The elastic energy is given by the product of the strain 
tensor Eij and the stress tensor Ty, integrated over the 
volume of interest: 
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Figure 1: A cylindrical mirror. 



which are Besscl's differential equations for r. The so- 
lutions are given by Bondu et al |15j with some typos 
corrected by Liu and Thorne [5]. Boundary conditions 
make a difference between the solutions with an infinite- 
size mirror and with the finite-size mirror. We will follow 
their calculation for the finite-size mirror and extend it 
with the coatings in Sec. [3] 



U = \ fj2 E v T v dV (h3 = r,^,z). (4) 2.3 Heat equation 



See Fig. [T] for the location parameters. The strain tensors 
of a cylinder with the axisymmetric pressure on the center 
are expressed by the displacement vectors u r and u z as 
follows: 

du r u r du z 

— , E, 



Either by adding the imaginary force or the imaginary 
entropy, the imaginary heat is generated and the finite- 
speed heat flow results in dissipation. The heat equation 
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and the stress tensors are as follows 
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Here A and \i are so-called Lamme coefficients: 
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with Y as Young's modulus and v as Poisson's ratio. The 
other elements of the tensors are zero due to the axisym- 
metry. The tensors should meet Newton's second law 
and Hook's law. Consequently, the elastic equation in 
the static case is summarized into two equations: 

9 V SV ldu r U r (8) 
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where 9j (z) is the time- varying temperature that is a 
Fourier transform of the temperature fluctuation 5T(z, t), 
Kj = fcj/Cj is the thermal diffusivity (fcj is the thermal 
conductivity and Cj is the specific heat per volume |17j). 
q]{z) is the heat source, and the subscript j indicates sub- 
strate (s) or coatings (c) . In the case of a coated material, 
the dissipation of the heat flow in the z-direction is larger 
than the r-direction, thus V in Eq. (fTT)l) shall be replaced 
by d/dz. The heat flow nCdO/dz gives the dissipation: 
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Here the bracket <> means the time average. 

The heat source q 3 is different between the equations 
for thermoelastic noise and thermorefractive noise. In 
the case of thermoelastic noise, the heat source is the 
expansion due to the imaginary force and it can be derived 
from the law of adiabatic temperature change [9] [18] : 
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where a- } is the thermal expansion, and Bj , the expansion, 
is expressed by the strain tensors: 
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In the previous studies [6]|10j. the expansion is regarded 
to be constant in z. This is good for the coatings while 
the z-dependence in the substrate, which is taken into 
account in our calculation, makes a non-trivial difference 
at low frequencies. We will explain the detail in Sec. |U 

In the case of thermorefractive noise, the heat source 
is the change of the refraction index caused by the imag- 
inary entropy perturbation [5] [5] : 
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where A is the wavelength of light, p(r) is the Gaussian 
profile of the beam, and /3 e ff is the effective temperature 
dependence of the refraction index given as follows: 



(16) 



with rii and n 2 as the refraction indices of two coating 
materials > n 2 ) and /?i and f3 2 as their temperature 
dependence. The delta function in Eq. (|15[) means that 
the heat source exists in the very beginning of the coating 
layers. See Ref. [5] for the details. 

As is shown in Ref. [6], the heat sources through the 
expansion and through the change of the refraction index 
have opposite signs in the heat equation. Both the ther- 
mal expansion a and the refraction-index change (3 are 
positive constants, i.e. the mirror expands geometrically 
and optically by the increase of the temperature, but the 
phase shifts due to the geometrical expansion and the op- 
tical expansion are opposite. Consequently, thermo-optic 
noise will be smaller than thermoelastic noise or ther- 
morefractive noise alone. We will show the calculation 
result with a finite-size mirror in Sec. [6] 

2.4 Mono-layer approximation 

The calculation of thermal noise in this paper, as well as 
other previous works, is based on the model that a single- 
layer coating with the thickness of multi-layer coatings is 
attached on a substrate. For Brownian thermal noise, the 
noise levels individually calculated with the silica coatings 
and with the tantala coatings should be square-summed. 
For thermoelastic noise, as is done by Fejer et al [10], we 
should replace some groups of parameters by the averaged 
value according to the following way: 
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where dj is the coating thickness of each material, d is the 
total thickness, the subscripts S and T indicate silica and 



tantala, respectively, and X, an operator to be averaged, 
would be the heat source g c , the thermal diffusion n c in 
Eq. ITU)) , or the thermal conductivity k c C c in Eq. (fT2"j) . 
For thermorefractive noise, /3 e ft in Eq. ([TFJ is already an 
averaged quantity. In this paper, we use a single tantala 
coating in Sec. H] and an averaged coating in Sec. [6] 

We should note that using a mono-layer coating is 
an approximation even with the averaging. Although the 
probe light is after all reflected by the coatings, some frac- 
tion of it transmits through the first few layers before re- 
flected by a later layer. More rigorous calculation would 
require solving the elastic equation and the heat equa- 
tion with as many boundary conditions as the number of 
layers, and it would probably include some coherent can- 
cellation of the volume fluctuation and the fluctuation of 
the refraction index. We shall leave this as a future work. 

In this paper, the optical length of each coating layer is 
a quarter of the wavelength of the probing light. Recently 
Principe et al has pointed out that thermal-noise level 
would decrease by tuning the layer thickness [2D]- We 
shall also leave this as a future work. 



3 Brownian thermal noise 

Let us follow Bondu's calculation to derive the strain and 
stress tensors in a cylindrical substrate. The tensors of 
coatings will be derived afterwards. The boundary con- 
ditions are 

T rz (r, z = 0) = 0, T zz (r,z = Q) = -F p{r) , 

T rr (r = a,z) = T rz (r = a, z) = , 

T zz {r,z = h)=T rz (r,z = h) = , (18) 



where 
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is the Gaussian profile of the beam. The solution to the 
elastic equation [Eqs. ©©] with the boundary condi- 
tions [Eq. HH1)] is 

u r = ^^A m (z)Ji{k m r) + Au r , (20) 

m 

u z = B m (z)J (k m r) + Au z , (21) 
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with 
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where Cm is a so-called Bessel-zero function that satisfies 
■MCm) = 0- and 
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The derivation is shown in Ref. [9] [15]. Without Am,. 
and Au z , the solutions (|20|) and (l2Tj) would satisfy all the 
boundary conditions but T rr (r — a, z) — 0. The addi- 
tional terms Au r and Au z make the difference approxi- 
mately zero. Actually these terms become dominant in 
the noise spectrum when we take the limit h <C a. 

Plugging Eqs. f20 |) (|2T ]l into Eqs. © © , we get the 
strain and stress tensors of the substrate. If we put them 
into Eq. (01 and then into Eq. © , substrate thermal noise 
would be calculated; the result is shown in Ref. [9]. 



With the coatings, as is introduced by Harry et al [5], 
the boundary conditions between the substrate and the 
coatings are: 
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(25) 



where the elements with a prime (') are for the coatings. 

(') (') 

At the boundary and in the coatings, ErJ and T r V are 
actually zero. Since the coatings are thin, we can assume 
that the strain and stress tensors are constant in terms of 
z. After some algebra, the strain tensor elements of the 
coatings are given as 
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then the stress tensor elements of the coatings are given 
as 



T' rr — (A' + 2fi')E' rr + X'{E'^ + E' zz ) , 
(X' + 2^)E^ + X'(E' ZZ + E' rr ) , 
= (A' + 2p!)E' zz + X'(E' rr + E'^) , 
T;, = 0. (27) 
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Putting these into 
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and then into Eq. ([3]), we obtain the power spectrum of 
coating thermal noise. Figure [2] shows the /i-dependencc 
and wo-dependence of the spectrum density \J S X (Q). One 
can see that the noise level agrees to the result with an 
infinite-size mirror, shown by dashed curves, with h suffi- 
ciently larger than ~ a, and increases by h~ 2 as the mir- 
ror becomes thin. The dotted curves in the top panel of 
Fig. [2] is the result with the thin-plate calculation, which 
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Figure 2: top: the /i-dependence, and bottom: the Wq- 
dependence of Brownian thermal noise. They agree to 
the previous results in the thin limit and the thick limit. 



we introduce in Appendix [A"l As h is sufficiently smaller 
than ~ a, the results with a finite-size mirror and with a 
thin plate coincide. Here the frequency is 100 Hz, the mir- 
ror radius is 2.5 cm, the beam radius in the top panel is 
1 cm, the thickness in the bottom panel is 2.5 cm, and the 
number of the coating layers is 3 for tantala and 2 for sil- 
ica; these are the parameters for a quantum-measurement 
experiment at Hannover |21j . 

Taking the limit a — > oo and h — ► oo, the strain and 
stress tensors agree to what are shown in Ref. [5] , and the 
noise spectrum is given as 
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of coating Brownian thermal noise for a gravitational- 
wave detector |23j and the result was identical to what 
we see in the bottom panel of Fig. [5J 

4 Thermoelastic noise 

The heat source of thermoelastic noise is the expansion. 
The expansion of the substrate is calculated to be 
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with B' m {z) — dB m /dz. The expansion of the coatings is 
calculated from Eq. (|26|) as 
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These expansions appear on the right side of the heat 
equation: 
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The homogeneous solution of the heat equation, which is 
the solution of the left-hand-side of Eq. (|32|) being zero, 
is 



9^ = Aj sinh (7jz) + B } cosh (7jz) , 



(33) 



where Aj and B s are coefficients that will be derived with 
the boundary conditions, and 7j is the complex propaga- 
tion constant given as 



7j 



(29) 




(34) 



For example, coating Brownian thermal noise of a mirror 
of a mirror in Advanced LIGO [22] , a second-generation 
gravitational- wave detector (a = 17 cm, ft, = 20 cm, wq = 
6.2 cm, and with 19 doublets of silica-tantala coatings), 
the thermal-noise level with the finite-size analysis is ~ 
2.6 % times smaller than that with the infinite-size anal- 
ysis. 

Coating Brownian thermal noise can be numerically 
obtained using the finite-clement method, although it takes 
longer time. Yamamoto et al calculated the wo-dependence 



The particular solution of the heat equation is given to 
cancel the right-hand-side of Eq. (|3"2")l . Reference [TU] uses 
an approximation that both S and C are constant in z, 
so that the particular solution is simply 

aP - a ' iY ' iT -6j (*~0), 



Cj(l-2uj) 

and the noise spectrum for the infinite-size mirror is de- 
scribed in an elegant form. In this paper, however, we 



6 



shall calculate the noise spectrum without this approxi- 
mation. The particular solution for the substrate without 
the approximation is 
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The particular solution for the coatings is simply 
oP _ a c Y c TF 



C c (l - 2v c ) 



(35) 



(36) 



with 9 C in Eq. (pTTj) . 

The complete solution is the sum of the homogeneous 
and particular solutions: 
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There are four boundary conditions to be met, which are 
(i) the heat flow at z — is zero, (ii) the heat flow at z = h 
is zero, (hi) the heat flows from the coatings to the sub- 
strate and from the substrate to the coatings are equal, 
and (iv) the temperature at the border is also equal: 



(i) k c C\ 



d9c 
dz 



= 0, 



(ll) K S C S — — 

oz 
09 

(hi) K C C C ^- 
oz 



= , 



z—d 



OZ 



z—d 



(iv) 



(38) 
(39) 
(40) 
(41) 



The coefficients in the homogeneous solution are then 
given as 



K = 



Ex/7. + (E 2 - 0P) 



cosh (7cd) + R sinh (j c d) 



where 



-(Ei/7 s )cosh( 7c d) + (E a - c p )i?sinh (j c d) d 



cosh ("fed) + R sinh ("fed) 



(42) 



R 



k c Cc7c 

K S C s 7s 



K S C S 



(43) 



and also Ei = 66^ /dz\ z= d ~ 86^ /dz\ z=0 and E 2 
Of (d) ~ C(0); namely, 



Ei 



a s Y s T 
C s (l-2u s ) 



E 



-iVik 2 n fi(a m +p m ) 
iCl — K s /c m A + 2/i 



Jo(k m r) 



~3A + 2^ ( 2Cl + 7i 



a s Y s T 
C s (l-2^ s 

1 

~3A + 2^i 



E 



% fa'/-, 



fj,(a m - /3 m ) 



ifl - K s /s m \ + 2fi 



(2c -po) 



Jo(k m r) 



(44) 



Here we use two approximations. One is to ignore the 
terms with e~ lsh in the presence of other terms in Eq. (|4"2"]) ; 
this is fine as far as the target frequency is higher than the 
inverse of the relaxation time of the temperature gradient. 
The other is to ignore the difference between e ±kmd and 
unity as well as to ignore the terms with d/h in Eq. (|44[) : 
this is fine if the beam radius is not as small as the order 
of d. Just in case, the followings are the terms that could 
be added to each term in Eq. (|4"2"]l : 



B c = — 



-7.(fc-«i) x [dOP/dz\ z=h ] 



7 S [cosh (7 C d) + Rsrnh (~f c d)] 
A s = _B c [sinh (j s d) cosh ("f c d) — -Rcosh ("f s d) sinh ("f c d)] , 
B s = 5 c [cosh (7 S d) cosh (j c d) — R sinh ("f s d) sinh (j c d)] . 

(45) 



Plugging the coefficients in Eq. (|42|) into the homoge- 
neous solution [Eq. ([33]) ]. adding the particular solution 
[Eqs. ([SSDdMD], and then putting 6 S and 6 C [Eq. (37])] into 
the following equation: 



W = 2tt 



K C C C 

2T 



-2tt- 



K S C S 

~2T 



dzrdr 



dz 



dzrdr , (46) 



we obtain the thermoelastic dissipating power, which with 
Eq. ([3]) gives the power spectrum of thermoelastic noise. 

Taking the limit a — > 00 and h — > 00, and ignoring 
the z-dependence of the expansion in the substrate, the 
noise spectrum agrees to the result of Ref . [10] : 



S x (Cl) 



16k B T 2 d(l + is s ) 2 a 2 s C c ~ 2 



nC 2 w 2 n 



A 2 (gi+g 2 ) 



(47) 



with 



~,Cs 



1 



1 



(1-2^ 



- 1 



2a s C c 1 — v c V 1 + v -- 
91 = (sinh£- sin 0/(0^) , 
<? 2 = i?(cosh£-cos£)/(£r D ) , 
T D = (l + # 2 )cosh£+ (l-# 2 )cos£ + 2iismh£ 
C = y/2nd 2 / Kc . 



(48) 



The noise spectrum with gi is the contribution of the 
dissipation in the coatings [the first term in Eq. (|4"6"]l ] and 
that with <7 2 is the contribution of the dissipation in the 
substrate [the second term in Eq. 
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z-direction, which is the more significant one for coating 
thermoelastic noise, we shall resolve Eq. (|4"9")l into two 
parts. Following the derivation of substrate thermoelas- 
tic noise shown in Ref. [9] , we can easily find out that the 
dissipation of the heat flow in the r-direction and that in 
the z-direction are equal, thus 



3 r 2 (l + v s ) 2 a 2 K; 



(50) 
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Figure 3: Thermoelastic noise at 100 Hz (top panel) and 
at 10 kHz (bottom panel) . The solid curve is thermoelas- 
tic noise of a finite-size mirror caused by the heat both in 
the coatings and the substrate. The black dotted curve 
is thermoelastic noise caused by the heat in the coatings 
and the gray dotted curve is the one by the heat in the 
substrate; both are with the thin-plate calculation. The 
dashed curves are thermoelastic noise of an infinite-size 
mirror. Substrate thermoelastic noise is calculated only 
from the dissipation of the heat flow in the z-direction. 



Figure[3]shows the /i-dependence of thermoelastic noise 
at two different frequencies. In the top panel, at 100 Hz, 
as h becomes small, the result with a finite-size mirror 
coincides to the level of substrate thermoelastic noise of a 
thin plate. We can see a non-trivial difference in a broad 
middle range between the result with our finite-size anal- 
ysis and previous results with the approximations. In the 
bottom panel, at 10 kHz, as h becomes small, the result 
with a finite-size mirror coincides to the level of coating 
thermoelastic noise of a thin plate. The result with a 
finite-size mirror coincides to the result with an infinite- 
size mirror as h becomes large. Substrate thermoelastic 
noise contributes more when the measurement frequency 
is low, the mirror is thin, and/or the coating is less. Note 
that the noise level with the finite-size analysis is larger 
than the square-sum of coating thermoelastic noise and 
substrate thermoelastic noise due to their correlation. 

5 Thermorefr active noise 

The heat source of thermorefractive noise is the change of 
refraction index. The heat equation for thermorefractive 
noise is 



A difference between our result and the result in Ref. [TU] 
according to the z-dependence of the expansion in the 
substrate could be compensated by adding substrate ther- 
moelastic noise of an infinite-size mirror derived with the 
thin-coating approximation in Ref. [3]: 



iil9 c 



8P 

' dz 



-in 



(3 cS XTF e" 



C c 



-p(r) , 



d 2 



(51) 



16k B T 2 {l + iy s ) 2 a 2 K s 

V^c s w*n 2 



(49) 



Here, we use the term substrate thermoelastic noise as 
what is caused by the heat source in the substrate due 
to the z-dependence of the expansion, while coating ther- 
moelastic noise is by the heat source at the interface of 
the substrate and the coatings due to the difference of the 
materials. Note that both are contained in the 6 S term 
and are dissipated in the substrate. Therefore, those two 
kinds of thermoelastic noise in the substrate should be co- 
herently summed up and our calculation offers a proper 
treatment by taking into account the z-dependence of the 
expansion in the substrate. Equation (|49[) is derived from 
the dissipation of the heat flow both in the r- and in- 
directions. Since our calculation contains only that in the 



where we have replaced the delta function to e~ z <^ with 
an arbitrary small length £, which does not appear in the 
final result as far as it is sufficiently small. The boundary 
conditions with a finite-size mirror are same as Eqs. 
(|4ip. The particular solution is 



-iCl (3 cS \TF e~ z/e 



in - k c /£ 2 
tnp cS ~\TF 



c c 



P(r) 



z/e p(r) 



. 



(52) 



and then the coefficients of the homogeneous solution, 
which is in the same form as in Eq. (|33p , are derived: 

= in p(r) , 

K c Cc7c 



8 



B r 



A, = iQ 



.^PcsXTFq 771 i? cosh (7 c d) 

k c Cc7 c 
/3cffATF 



772 sinh (7 c rf) and the particular solution for the substrate is same as 



k c C c 7c 



p(r) 



where 



— U 



?7o = 
»7i = 

»/2 = 



771 i? sinh (j c d) + r\i cosh ("fed) 'what we have derived for thermoelastic noise [Eq. ([55]) ]. 

^ The boundary conditions with a finite-size mirror are 

same as Eqs. Let us use the approximations 

(|45| and (|54|) . which have been proven to be safe in the 
(53)individual calculations for thermoelastic noise and ther- 
morefractive noise. The coefficients of the homogeneous 
solution are then given as 



772 cosh (7 C d) + rji R sinh (7 C d) 



I - cxp (-2 7s fe) 
1 + cxp (-2-f s h) l ~ J 
cosh (7 S d) — r/o sinh (j s d) 
770 cosh (7 S d) — sinh (7 s <i) 



(54) 



As 770 — f, thus 771 ~ 772, Eq. (|53[) can be more simplified. 
In fact, after the simplification, none of the coefficients 
contains h. The thermorefractive-noise level with a finite- 
size mirror is then almost same as that with an infinite- 
size mirror: 



S x (fi) = 



2\72fc B T 2 r 2 I I 



K c Cc TTWq 



/3 c 2 ffA 2 



(55) 



where 



T 2 = (I + R 2 ) sinh£ + (I - i? 2 ) sin£ + 2i?cosh£ , (56) 

and Yd has been given. We introduce earlier than To 
and Ti, which will be shown in Sec in order to keep 
the same notation as Ref. [BJ. 



6 Ther mo-optic noise 

It has been pointed out by Evans et al that thermoelas- 
tic noise and thermorefractive noise should be coherently 
added with a proper treatment due to their common ori- 
gin; it is now called thermo-optic noise [BJ. In Sec. H] and 
Sec. [5j we have studied the behavior of these two kinds 
of noise and derived the individual noise levels, but what 
should be used for the noise estimation is the result in 
this section. The heat equation of thermo-optic noise has 
both thermoelastic heat source and thermorefractive heat 
source on the right-hand side: 



ine c 



= -in 



92 



p(r) — in 



c c t 

z 



a c Y c T 
C c (l - 2v c ) 



= -in 



C s (l-2v s ) 



6s 



(57) 



The homogeneous solution is in the same form as Eq. ([55)) . 
The particular solution for the coatings is 



inf3 cB \TF a _ ; 
e 



It 



p(r) - 



a c Y c T 
C c (l - 2i> c 



(58) 



B c 

A ~ - 

where 

ni 



n. 



ill — p(r) , 



K cC c 7 c 

S1/7. 4 



(£2 + Hi) 



cosh ("f c d) + R sinh (7c d) 
-(£i/7 s )cosh ( 7c d) + (Ej, 



cosh ("f c d) + i?sinh (7 c d) 



IEJ-Rsinh (7 C d) d 



-A 



(59) 



OL c Y c T 



cr c (i-2i/ Lj 

ini3 cS XTF 

K cCc7c 

a c Y c T 
C c (l-2u c ) 



p(r) [sinh (7 c d) + R cosh ("f c d)] 
inp eS XTF Q 



e 



n c C c "f c sinh (7 c d) 



p(r) . (60) 



One can see that the coefficients in Eq. (|59"|) coincide to 
those in Eq. (|42| if a thermorefractive constant /3 e ff is 
supposedly erased and they coincide to those in Eq. (|53|) 
if thermoelastic constants a c and a s are erased. 

As well as we have done in the previous sections, plug- 



ging 8 S and 9 C with the coefficients in Eq. 
following equation: 



into the 



W = 2tt 



K C C C 

2T 



-2tt 



K S C S 

2T 



dzrdr 
' d6 R 



dz 



dzrdr 



we obtain the thermo-optic dissipating power, which with 
Eq. (PJ gives the power spectrum of thermoelastic noise. 
Taking the limit a — > 00 and h — > 00, and ignoring the 
z-dependence of the expansion in the substrate, the noise 
spectrum agrees to the result of Ref. [6J : 



s x (n) = 



2k B T 2 1 I I I 
n ^fCc^rif^ 



r (Aad) 2 - TxAad • /3 eff A£ + T 2 (f3 cff Xtf 



where 



Aa = 2a s (l 



-4 A 



(61) 



(62) 



9 



and 



r = 2(sinh£-sin£) + 2i?(cosh£-cos£) 



f f £ 
! sin - R cosh - + sinh 
2 I 2 2 



r 



See Eq. (@SJ) for A, and Eqs 
respectively. 

With the averaging, some of the constants in Eq. (JBT 
should be replaced as follows: 



(63) 

SBJ for T D and T 2 , 



Aa — > Aa = a c — 2a 8 (l + v a ) 



a 



with 



and 



1 + v s fl + v, 



1 — V r V 1 + V. 



(l-2*)£ 



(64) 



(65) 



avg 



C c = Cs -f + C T 



i 

cv 



1 d s 



i d T 



ksCs d ktCt d 



(66) 



then £, 7 C , and i? should be replaced to the averaged 
constants with these new C c and k c . 

The averaging can be done for the calculation with a 
finite-size mirror as well. With a proper treatment, III 
and II2 are replaced to 



1 - 2z/ ( 

inp cS \TF 



p(r) [sinh (*y c d) + R cosh (7 c d)] , 



1 



1 - 2i/ c 



KcC , c 7 c sinh (7 c d) 



p(r) . 



(67) 



Replacing C c , k c , £, 7 C , and R to the averaged constants 
shown above as well, we obtain the noise spectrum with 
the multi-layer coatings of adequate accuracy. 

Figured] shows the noise spectra of thermo-optic noise 
with a finite-size mirror and an infinite-size mirror. Here 
the averaged coatings are used. Compared with the re- 
sult with an infinite-size mirror [6], thermo-optic noise 
with a finite-size mirror is larger at low frequencies and 
smaller at high frequencies according to the difference in 
thermoelastic noise, while thermorefractive noise is same. 

It is worth noting that there is a small, conceptual dif- 
ference in the treatment of the two kinds of heat source for 
thermo-optic noise between this paper and Ref. [6|. Rig- 
orously speaking, it is imaginary force that should be in- 
jected to calculate thermoelastic noise and it is imaginary 
entropy to calculate thermorefractive noise. In Ref. [B] it 
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Figure 4: Thermo-optic-noise spectrum of each mirror 
in a Fabry-Perot arm cavity of Advanced LIGO detec- 
tor [55] • The mirror radius is 17 cm and the mirror thick- 
ness is 20 cm. The beam radius is 5.5 cm on the input 
test mass and 6.2 cm on the output test mass. The input 
test mass has 8 doublets, and the end test mass has 19 
doublets of silica-tantala coatings. 



may seem like the entropy is injected for both purposes, 
and in our paper it may seem like the force is injected for 
both purposes. In fact, which imaginary quantity is in- 
jected does not matter since the conversion from the tem- 
perature fluctuation to displacement noise on the surface 
can be simply done by multiplying some constants. This 
is true as far as the stress inside the coating is uniform, 
or in other words, all the light is regarded to be reflected 
at the surface of the first coating layer. 

7 Summary 

We demonstrated the calculation of Brownian thermal 
noise and thermo-optic noise (thermoelastic noise + ther- 
morefractive noise) of the coatings on a finite-size cylin- 
drical mirror. We used a method developed by Bondu et 
al based on the fluctuation-dissipation theorem to calcu- 
late the elastic response of the mirror, and extended it 
with the coatings. Comparisons with the previous calcu- 
lations of thermal noise with an infinite-size mirror and 
with the independent calculation with a thin plate val- 
idate our results with a finite-size mirror. We showed 
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how Brownian thermal noise and thermoelastic noise in- 
crease if the mirror becomes thin. We treated all the heat 
sources in the coatings and the substrate coherently, and 
the results give us the most accurate estimate of coating 
thermal noise. 
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A Thermal noise of a thin plate 

We calculate thermal noise of a thin plate based on the 
study by Yamamoto et al [24] . The noise level can be cal- 
culated with the modal-expansion method easily as the 
contribution of higher order modes is extremely small in 
the case of a thin mirror. This method is quite indepen- 
dent from the method we have shown in the main body 
of this paper, and the agreement of the results validates 
the calculation with a finite-size mirror. 

A.l Brownian thermal noise 

With the modal-expansion method, the noise spectrum 
of the thermal motion in the first mode is given by 



s x (n) 



4k B T 1 



(68) 



where mi and uj\ are the effective mass and the resonant 
frequency of the first mode, respectively, and Q c s is the 
effective Q value of the coatings, which can be derived 
from the intrinsic loss angle and a compensation factor 
to adjust the elastic energy in the coatings to the total 
elastic energy [2"5] : 



1 

QcS 



3Y c d 
~Yh' 



(69) 



In the case of the thin-mirror analysis, contributions from 
the higher order modes are negligible, so that we can just 
calculate the contribution from the first mode. The elastic 
equation is 



h 2 Y R 



12(1 - v») 



A 2 n7i(r) = -p s w\wx{r) , 



(70) 



with tJ7i(r) as the one-dimensional displacement of the 
first mode and p s as the density of the substrate. The 



boundary conditions are 
d (d 2 



1 dw\ 



dr V dr 2 r 
l w\ v s dw\ 



dr 



+ 



dr 2 r dr 
which gives the resonant frequency as 







Y s h 2 



12A(1-". 2 ) 



(71) 



(72) 



where ot\ = 2.9493 is a solution of 
2(1 -V,) Tt ^ Tl ^ Jl{*l) 



Oil 



-Ji(ai) - Jo(ai) 



Ii(ai) 



Jo(ai) = , (73) 



with J n {a\) and I n (oi) as the Bessel function and the 
modified Bessel function of the first kind, and then z&i is 
given as 



j ( n -Mgi) j ( 

w± = J I ai- I - — — r-io a 

V a/ hyoii) V 



The effective mass is calculated from 



(74) 



mi 



J p s \mi\ dV 



Wip(r)dS 



(75) 



Plugging Eqs. 1[59" 1) (|72" 1) ([75 ] ) into Eq. |B5]), we obtain the 
spectrum of Brownian thermal noise in the coatings of a 
thin plate. 

A. 2 Thermoelastic noise 

Let us first derive thermoelastic noise caused by the heat 
in the coatings. Here we use the thin-coating approxima- 
tion, and also we assume a mono-layer tantala coating. 
The solution of the heat equation can be resolved into 
the sum of functions that meet the boundary conditions 
that the heat flow should be zero at z = and z = h, 
then, 



=5> 



2 / nwz \ 
— cos I I 



h~~~\ h J ' 

n 

Plugging this into the heat equation, we get 

a e sY s Td 



(76) 



a c ff 
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ifl + K s (nir/h) 2 C s (l - 2v s ) 
Y c {\ - v B ) Cc 

a °Y s (l-v c ) as C s ' 
1 — 2v s h / d 2 zui 1 dwi 
1 — ^ s 2 V dr 2 r dr 



9 x 



h 1 



(77) 
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The dissipation power is then given as 



W = 



K S C S 

2T 



89* 



dV , 



(78) 



and the total energy is 



24(l + i/.)(l 



dr 2 
2(1- 



- Vs) 

1 dw\ 
r dr 



rdr 



d 2 vj\ 1 dw\ 
dr 2 r dr 



rdr 



(79) 



As the lost energy in one period is the total energy mul- 
tiplied by 2%(p, the loss angle of the first mode caused by 
the heat source in the coatings is given as 



i coa 
°1 



w 



uxEl 01 
Y s a 2 s T 1 



6d 2 



C s 

where r = h 2 /(k s it 2 
d 2 w\ 



1 - 
and 



h 2 



fir(n) 



(fir) 2 + (n)' 



(80) 



B x = 



1 dw\ 



dr 2 



r dr 



rdr 



d 2 w\ 



dr 2 



1 dw\ 
r dr 



rdr 



-2(1 - i/O / 



d TJ7l 1 cfol 



dr 2 r dr 



rdr 



,(81) 



which is numerically calculated to be 1.47232. Substitut- 
ing l/Qcff in Eq. (|69j) to </>i° a , we obtain coating ther- 
moelastic noise. Note that the thin-mirror approxima- 
tion let some errors in the result with a thin plate at 
high frequencies due to the thin-coating approximation; 
for example at frequencies higher than ~ 1 kHz with 
a = 17 cm, h = 5 cm, w$ — 6.2 cm, and N — 19. It 
is not a problem in Fig. [3] as the number of coatings is 
only 3. 

For thermoelastic noise caused by the heat source in 
the substrate, the loss angle of the first mode is given as 



6 »ub 



Ym 2 T 1 



C s l-i/ B l + (fir) 



(82) 



the derivation of which is shown in Ref. [27j . Substituting 
1/Qcff in Eq. (f69|) to </>f lb , we obtain substrate thermoe- 
lastic noise. 



B List of the parameters 

• Coating loss angle 0j 

Si0 2 : 1.0 x 1CT 4 , Ta 2 5 : 4.0 x 1CT 4 



• Thermal conductivity fcj 

SiOz : 1.38 W/m-K, Ta 2 5 : 33 W/m-K 

• Thermal expansion ay 

Si0 2 : 5.1 x 1CT 7 /K, Ta 2 5 : 3.6 x 1CT 6 /K 

• Specific heat per volume Cj 

Si0 2 : 1.64 xlO 6 J/K-m 3 , Ta 2 5 : 2.1 xlO 6 J/K-m 3 

• Thermal diffusivity kj {— k-JC- s ) 

• Young's modulus Yj 

Si0 2 : 7.2 x 10 10 N/m 2 , Ta 2 5 : 1.4 x 10 11 N/m 2 

• Poisson ratio Vj 

Si0 2 : 0.17, Ta 2 5 : 0.23 

• Refraction index nj 

Si0 2 : 1.45, Ta 2 5 : 2.06 

• Temperature dependence of the refraction index (3j 
Si0 2 : 8 x 10- 6 /K, Ta 2 5 : 14 x 10" 6 /K 

• Density p 

Si0 2 : 2200 kg/m 3 

• Wavelength of light A 
1064 nm (Nd:YAG laser) 

• Temperature T 
300 K 
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